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ABSTRACT 

It was suggested by several authors that hypothetical primordial black holes (PBHs) may 
contribute to the dark matter in our Galaxy. There are strong constraints based on the 
Hawking evaporation that practically exclude PBHs with masses m p bh ~ 10 15 — 10 16 g 
and smaller as significant contributors to the Galactic dark matter. Similarly, PBHs 
with masses greater than about 10 26 g are practically excluded by the gravitational 
lensing observation. The mass range between 10 16 g < m p f,h < 10 26 g is unconstrained. 
In this paper, we examine possible observational signatures in the unexplored mass 
range, investigating hypothetical collisions of PBHs with main sequence stars, red giants, 
white dwarfs, and neutron stars in our Galaxy. This has previously been discussed as 
possibly leading to an observable photon eruption due to shock production during the 
encounter. We find that such collisions are either too rare to be observed (if the PBH 
masses are typically larger than about 10 20 g), or produce too little power to be detected 
(if the masses are smaller than about 10 20 g). 
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Introduction 



The idea of dark matter (DM) in the Universe was first advanced in the 1930s, when lZwickvl (|1933l . 
19371 ) discussed the deviation of the masses of galaxy cl usters away from expected values. The ide a 
was subsequently confirmed in many different ways (e.g. lOstriker et al.lll974l : lKomatsu et al.ll2009i ). 
The strong deviation of the rotation curves of spiral galaxies away from what would result from 
just the luminous matter, has given an important addition al quantitative indication o f the existence 
of DM. Its presence in our Galaxy was first discusse d by iKahn and Woltjer The current 

estima te gives Mum = 9 • 10 11 M Q ( Xue et al.1 l2008l ). which matches that by IKahn and Woltjer 
(|1959l ). Searches for non- luminous DM in the form of brown and white dwarfs, neutron stars and 
black holes, with masses above 0.1 — 0.2- Mq ("MA CHOS") has shown that these objec ts are far too 
rare to make up a significant part of the total DM (jAfonso et al.ll2003l ; lYoo et al.ll2004l ) . Within the 
past five years, it has even been shown by WMAP that baryonic matt er only makes up ~ 4% , of the 
critical density of the Universe while dark matter makes up ~ 20% (IKomatsu et al.ll2009i ). While 
objects like neutron stars or white dwarfs are thereby excluded as making up a major proportion 
of the DM, black holes with masses between 10 16 g < m p bh < 10 26 g might still be responsible for 
a significant fraction, as explained in the next section. In certain models, these primordial black 
holes (PBHs) are predicted to be produced at an interesting level in early phases of the Universe 
due to density fluctuations of the type also observed in the cosmic microwave background. Their 



gravitational radius r g would be subatomic, 



10~ 8 • (m pbh /(10 20 g)) cm. The Bondi-Hoyle 



accretion radius, on the other hand, can be significantly larger. 



In this paper, we check whether hypothetical PBH collisions with Galactic objects may be ob- 
servable. A detection of observational signatures of such collisions would support the hypothe- 
sis that PBHs may represe nt a sign if icant fraction of the DM in the Galaxy. We find that the 
PBH collisions discussed by IZhilvaevI (|2007l ) would produce no presently-observable signatures. To 
prove our case, we show that even with the most optimistic conditions assumed, hypothetical col- 
lisions of PBHs with Galactic objects place no observational constraints in the PBH mass range 
~ 10 16 g < mpbh < 10 26 g considered in this paper. 

The assumptions made are as follows: 



• The mass distribution of PBHs is a delta-function, i.e. all PBHs have the same mass. The mass 
is treated as a free parameter, and the whole range ~ 10 16 g < m p bh < 10 26 g is considered. 
Any other distribution would significantly decrease the event rate: if we assume that all PBHs 
have the same mass and then consider this mass as a free parameter, we maximize the event 
rate for each single mass value, as we assume that 100% of the total mass is focused at one 
arbitrary value. For wider mass distributions, the total mass of the PBHs is less than 100% 
for a given mass value, which always reduces the event rate. 



• The entire energy from a single interaction is taken to be in the energy band of the potential 
detector. In reality, we expect the signal to be spread from radio wavelengths up to soft 
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gamma-rays. Using the entire radiated energy therefore gives an upper limit for detectability. 
A single photon is sufficient for detecting the event. Even for a b ackground-free sour ce, three 



photons are required for giving 90% confidence significance (see (jAmsler et al.ll2008l )) and so 
with a more realistic estimate, the detection probability would be further reduced by a factor 
of between 3 and more than ten. 

The signal is radiated at full opening angle, i.e. 2tt. A more realistic estimate would include 
consideration of the shock opening angle, reducing the expected observation probability. 

We assume that the signal is always strong enough that it can reach us over the entire thickness 
of the Galactic disk, h f» 300 pc. 



Further, we make a few assumptions that greatly simplifies calculations but could only slightly 
(indeed insignificantly) influence the final result: 



• Luminous matter is taken to be distributed uniformly in the Galactic disk and dark matter 
is taken to have constant density in the spherical halo. We will show that the enhanced 
density in the Galactic Center yields slightly better results, but with the same conclusion 
that detection of the events can be neglected. 

• We use a delta-function for the distribution of relative velocity between the luminous and 
dark matter, vo = 2.2 • 10 cm/s. We will show in this paper that the reduced velocity at the 
central part of the Galaxy does not influence our result. As for the possibility of a capture at 
low velocities, we show in Section [5] that this does not occur within the energy range being 
considered. 

• The PBH is taken to be uncharged. 



For alternative approaches to explaini ng dark matter by extending the standard model of particle 

physic s to include new particles, see e.g. lHaber and Kanei ()1985l ) ; iBertone et al.l (|2005l ) ; iBiermann and Kusenko 
(|200fih : lHooper and Profumcl d200^ . 



In order to investigate the interaction of PBHs with objects in the Galaxy and possible observational 
consequences, we need to ask the following questions: 



1. How frequent are the events? 

2. How energetic is each single event? 

3. What are the possible signatures? 



In the following sections, we answer each question separately for the five most important object 
classes in the Galaxy, i.e. main sequence stars, red giants (in particular their dense cores), white 
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dwarfs and neutron stars. The main intrinsic properties of those objects are listed in tableQ]together 
with their abundances in the Galaxy. Both the interaction probability and the energy loss strongly 
depend on the properties listed. Also, we consider collisions with the Earth, in order to investigate 



the cl aim that the Tunguska event might have been caused in this way (see iJackson and Ryan 



41973)). 

This paper is organized as follows. Section [2] summarizes the already-existing observational con- 
straints. In Section [31 interaction probabilities in the Galaxy and energy loss of a single event 
are derived. Possible signatures and the detectability of such events are investigated in Section [H 
Section [5] discusses the possibility of capturing the PBH in a stellar object. Finally, Section [6] 
summarizes the conclusions to be drawn from our calculations and discusses further possibilities 
for detection requiring additional assumptions. 



Primordial black holes - constraints 



Ever since the hypothesis of primordial black holes was first established ([Hawking] Il97ll ) . different 



theoretical and observational constraints on their actual abundance in the Universe have been 



discus s ed. While theore tical constraints are rather speculative and model dependent (see e.g. ICarr 
(|2005l ): iKhlopovl (120081 )). the observational ones are far more robust. They firmly constrain the 
abundance of primordial black holes in two mass ranges. In this section, we briefly review them. 
Since we are interested in PBHs as a possible contributor to the dark matter, we only consider 
PBHs that can have survived until today, that is those with initial mass above 5 • 10 g. 

There are three main methods of looking for indirect detection of primordial black holes above 
5 • 10 14 g. The first one aims for the detection of PBH evaporation. This method is sensitive to 
the lower mass limit of the initial mass function (IMF). This initial mass function is expressed 
in terms of the number of PBHs at the ir initial mass m, as a function of the initial mas s itself, 
dn/dm. As predicted by iHawkingi (|1974l ). and shortly afterwards also by iBekensteinl (|1975l ). black 
holes evaporate due to quantum fluctuations around the Schwarzschild radius of the black hole, 
leading to the evaporation of PBHs with masses below 5 • 10 14 g before the present time. Therefore, 
non-detection of the high-energy photons which would be produced by black hole evaporation, i.e. 

> 20 MeV, leads to a robust constraint on the mass function around 10 15 g. At 10 15 g, EGRET 
observations between 30 MeV and 120 MeV imply t hat the PBH mass density must be smalle r 
than 3.3 • 10~ 9 times that of the total dark matter (jSreekumar et al.l Il998l ; iBarrau et al.ll2003l ). 
Between 10 15 — 10 16 g, constraints from EGRET measurements are still at a level of < 10 -2 and 
provide significant constraints as well. At si gnificantly higher m as s scales, micr o - lensi ng can help 
to restrict the abundance of PBHs, see e.g. lAfonso et al.l (|2003l ). lAlcock et al.l (|1998l ) looked for 
planetary mass dark matter via gravitational micro-lensing in the direction of the Large Magellanic 
Cloud, resulting in constraints of up to 10% in the mass range 10 26 g < m pbh < 10 33 g. At even 
higher mass scales, 10 33 g < m p bh < 10 40 g, results from anisotropy and spectra measurements of 
the Comic Microwave Background (WMAP-3years and COBE-FIRAS) can be used to constrain 
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the number of PBHs (Ricotti et al 



200a). The excluded regions for the mass function are shown 



in Fig. [TJ In addition, iGouldl (|1992l ) proposed using gamma ray bursts (GRBs) to detect "femto- 
lenses" produced by dark matter objects in the mass range of 10 17 g < m p bh < 10 20 g. This method 
has not yet been applied, however, and so no experimental limit has yet been set. A mass range of 
about ten orders of magnitude, 10 15 g < m p bh < 10 26 g, then remains at present unexplored. 

In the standard scenario of density p erturbations in the early Universe, the IMF of PBHs, dn/dm(m), 
behaves as, e.g. lHalzen et al.l (|199ll ) 



dn 
dm 



oc m 



-2.5 



(1) 



and so most PBHs have masses at around the current evaporation level. The constraint by EGRET 
at the lower end of the mass scale therefore has significant consequences for the abundance of black 
holes at higher masses. 



The above model for PBH production is not, however, the only one. Several other proposed 
theories for PBH production predict significant contributions from PBHs at different mass scales. 
Examples are the collapse of topological defects, softening of the equation of state and q u antum 
gravity s c enari os, which, although speculative, deserve to be tested (see reviews bv ICarrJ (|2005l ): 



Khlopovl (|2008l ) and references therein). 



3. A toy model to estimates of the collision probability and energetics 

In order to describe possible observational signatures of the hypothetical collisions of PBHs with 
stars in our Galaxy, one must answer the two fundamental questions: how frequent are the colli- 
sions? and how energetic is a single collision? Precise answers to these questions are obviously very 
difficult to obtain, as they depend on genuine astrophysical uncertainties that cannot be resolved 
today. In the view of these uncertainties, we estimate the collision frequency and energetics using 
a very general and rather simple analytic toy model that adopts several simplifying assumptions. 
For example, we assume that all PBHs have the same mass m p bh- We do not specify that value. 
Instead, we consider m p bh to be one of a few free parameters of the model. 

Adopting ideas from the kinetic theory of gases, one may argue that a reasonable estimate of the 
collision frequency between stars and PBHs should depend only on the size r of the region of 
interest (e.g. either the whole Galaxy, or its central bulge), numbers of stars N* and PBHs N p bh in 
this region, the relative velocity v* between stars and PBHs, and the collision cross section o~ p bh,*- 
Similarly, a reasonable estimate of the collision energetics should depend only on the mass and 
radius of the star, M+, R*, mass of the black hole m p bh, and the relative velocity 

The final results obtained this way fully justify the adopted strategy. They show that the hypo- 
thetical collisions would be far too infrequent, or far too energetically weak, to be detected. This 
holds in the whole range of m p bh and all other parameters of the model. The mismatch involves 
several orders of magnitude. Thus, our conclusion that there are no observational constraints from 
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hypothetical collisions of hypothetical dark halo primordial black holes with Galactic objects cannot 
be changed by including fine details, for example by considering some particular mass functions for 
m P bh- 



3.1. How frequent are the collisions? 

Let us consider a spherical part of the Galaxy with the radius r containing N±(r) stars. We will 
consider stars of a specific type, e.g. main sequence, white dwarf, or neutron star. The dark matter 
mass in the region is MoM{ r )- The number density of primordial black holes can be estimated by 
assuming that PBHs represent a fraction r/ of the dark matter in the region, and that the mass of 
a single PBH is m p bh (v an d m p bh are free parameters of the model), 

r}-M DM (r) 1 

n p bh{r) = . 5 • (2) 

47T/3H rripbh 

The collision frequency is given by, 

n(r) « iV*(r) • j pbh (r) ■ a pbh ,*(r) « N*(r) ■ V 'Jfj 3 ^ ■ ^ • a pbh ^(r) . (3) 

4/3 7T H m pb/l 

Here o~ p bh,+ is the cross section of the interaction, and j p bh( r ) = n pbh( r ) ■ V( r ) is the flux of PBHs 
that move with respect to stars with the relative velocity V(r). We approximate V(r) by the 
rotational velocity in the Galaxy, which is close to radial. 

The cross section o~ p bh,*(r) = ir R${r) should be calculated from the energy and momentum con- 
servation during the encounter of gravitationally interacting star and PBH. The impact parameter 
Ro{r) of the encounter is defined by the PBH angular momentum, equal to 

m p b h -V(r)-Ro(r). (4) 

As we are interested in an optimistic estimate of the event rate, and hence of the cross section, 
here we use the conservation of angular momentum in order to put an upper limit on the impact 
parameter Rq for which there is a collision. This means that we use the maximum value for Rq 
for each collision, even if this value may be smaller in reality. This maximum value for the impact 
parameter corresponds to the limiting case of a PBH only grazing the star's surface, so that from 
angular momentum conservation we get 

m pbh V(r)Ro = m pbh v+R± R = R *y^ ( 5 ) 

where v± is the final PBH speed, that can be computed by considering energy conservation in the 
limit of m p bh <C M*: 

1 2 GAh 1 2 GM, 

2 V "7i^ = r*-^r (6) 

with M* being the mass of the star and < R > being the typical distance between stars. The char- 
acteristic distance at which the PBH enters the gravitational field of the star can be approximated 
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as the typical distance between objects in the Galaxy, assuming that each star-like object deflects 
the PBH at least marginally. We take this typical distance to be < R >~ 1 pc. We find that the 
gravitational energy is negligible in the initial stated, so that the final velocity can be expressed as 

/ o C M 

v, « sjv\r) + = vV 2 (r) + «Lc, (7) 

where v esc = 2GM t /^ is the escape velocity of the star. Combining equations (|T|) and (0) one 
arrives at 

p ( r \ — d v *( r ) _ R V^ 2 (r)+^ c 

^ o(r) " ^ ■ 7(0 " W) ' ( } 

Replacing the impact parameter Rq in Equ. ([3]) with the help of Equ. ([8]) we get, finally 



4/3 r 3 m pbh V^( r ). 

Now, the radial dependence of the parameters needs to be specified. The observed velocity profile 
of the Galaxy can roughly be approximated as V(r) = vq ■ r/rjj for < r < ~ 2.5 kpc, and 
V(r) = vq for r<2 < r < r^aio ~ 50 kpc. Starting from this observational fact, mass densities and 
masses enclosed in a radius r can be determined, as listed in table [21 This naturally divides our 
estimates in a central part, considering only the bulge of the Galaxy, and an outer part, considering 
the part of the Galaxy where the velocity is constant. Here, we will first consider the outer part of 
the Galaxy. In Section [3.31 we wm show how our results change for the central part. 

If the region considered is the whole outer part of the Galaxy (the "flat rotational curve region" 
outside the bulge), one may use r = rh a i ~ 50 kpc, M^m ~ 9 • 1O 11 M , and V(r) = vq « 
2.2 • 10 7 cms -1 . Introducing units convenient in this region, JV*g = N+/W 9 , t^o = v */ v o, R*,9 = 
-R*/(10 9 cm) and m p bh,20 = "V>/i/(10 20 g), one writes ([9]) in the form 

-3.„-l at _ „ „,2 D 2 f V \ _-l 

2.2 ■ 10 7 cm/ S/ 



h « 2 • lO^yr" 1 • JV*,9 ■ V ■ <o ■ ^,9 • ( 7717^7377: ) " ^1,20 ■ ( 10 ) 



Apart from the PBH mass (which is a free parameter) this estimate depends only on the stellar 
properties. For instance, white dwarfs (with v+/vq = 51/2.2) have an event rate of 1.4 yr _1 for a 
black hole mass of 10 20 g. 



3.2. How energetic is a single event? 

In estimating the kinetic energy which is transferred from the black hole to the star and then 
radiated, we follow previous studies of the motion of accretors in a continuous medium, see 



1 For < R >— 1 pc, we get GM+j < R >~ 7 ■ 10 3 cms -1 . The initial gravitational potential only becomes 
important for distances below < R >< 0.01 pc, where the term GM+j < R > becomes comparable to the initial 
velocity. Since such short distances between stars are rather rare, we neglect this effect. As a comparison, the distance 
between the Sun and its closest neighbor is ~ 1.3 pc. 
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e.g. iRuderman and Spiegel! (|197ll ). In making this first estimate, we assume a uniform and fric- 
tionless medium and disregard self- gravity of the gas. 

The black hole passing through the gas focuses matter behind itself, creating a wake, because of 
the gravitational interaction transferring momentum to the atoms nearby. The inhomogeneity in 
the medium gives rise to a force acting on the black hole, causing a decrease in kinetic energy. 



As pointed out by lOstrikerl (|1999l ) , the force on the black hole (and hence the energy loss) depends 
on whether the motion is supersonic or not, as well as on the degree of supersonic motion. In 
making a first order approximation, we consider the black hole to be moving at constant speed. 
Also, we assume the velo city to be supersonic, so that the drag force, and hence the energy loss, 
reduces to that given by IRuderman and Spiegell (|197ll ). The energy loss due to the bow shock is 
negligible compared to that due to the tail shock, and the energy loss per time is then given by 



dE loss ^G 2 m 2 bh p+ D 



dt 



hi 



D n 



(11) 



where D mSuX is the maximum linear dimension of the star (i.e. its diameter) and .D mm is the linear 
dimension of the accretor, lying between the black hole's horizon radius and its Bondi-Hoyle accre- 
tion radius. The exact formulation of the linear scales is not crucial for this first estimate, as they 
appear in a logarithm. Here, we use ln(Z) max /.D mm ) = 10. The total energy loss is then given by 



Eh 



dE lc 



loss 



dt 



dt 



In 



(12) 



considering a time interval of dt « 2 ■ R+fv*. Equation (fT2|) can be expressed as 

-2 



E\ 



2 ■ 10 27 erg • m 



'°' R *' 9 ' \ 2.2-107 C m/s 



• P*,5 



(13) 



where p* t s = /^/(lO 5 g cm 3 ). Table [3] lists the energy released and the interaction probability for 
black holes of 10 20 g passing through the different types of Galactic objects. 

In our picture, the shock front will locally heat the plasma to an extent depending on the particular 
Galactic object under consideration. The spectrum is therefore a thermal one. In the following, we 
discuss the peak energy of the spectrum in order to estimate the detectability for detector systems 
observing in the relevant energy range. 

To calculate the temperature of the shock, we assume supersonic shock conditions and a perfect 
gas, giving 

T ^^ Hi! >^-^i^w) ■ (14> 

where the downstream ve locity is diminished by a f a ctor of four with respect to the upstream 
velocity (see, for example. lLandau and Lifshita (119871 ): iParksl ()199ll )). Here, ks is the Boltzmann 
constant, is the velocity of the shock, run is the mass of a hydrogen atom, M is the Mach 



- 9 - 



number, and the average relative velocity between primordial black holes and stars is taken to be 
2.2 • 10 cm s _1 . Wien's law now gives the frequency at maximum intensity as 

£5 = £ • l^S-vSM* > 4.5 • 10- Hz • ( — -S A 2 , (15) 



6 b 64 jfc B V 2 - 2 ' 10 7 cm/s / 

where b = 2.9 • 10~ 3 Km. This can also be expressed in terms of the peak photon energy 

^neak = hVrnax > 0.2 keV • ( % -] . (16) 

p V2.2-10 7 cm/sy 

This is merely a lower limit, but should suffice as a first-order estimate. For a collision with Earth, 
an extreme ultraviolet spectrum is expected, while collisions with main sequence stars and red 
giant cores are expected to give an X-ray spectrum. White dwarfs and neutron stars cannot be 
treated with simple thermodynamics: we here assume that the shock velocity v* determines the 
temperature of the shock. Thus, we get a peak energy which is increased by a factor 4 2 = 16, 

£neak = 1-6 keV ■ ( % -] . (17) 

p \ 2.2- 10 7 cm/ S y v ; 

For a white dwarf, we get a soft 7 spectrum, while for a neutron star we get high-energy gamma 
rays, see table El Since we are dealing with order-of-magnitude estimates here, the actual peak 
energy could deviate from what is calculated. Therefore, we present our results for all objects and 
for all possibly interesting energy ranges, i.g. X-rays up to 7— rays. 



3.3. A scaling to the bulge 

The estimates of frequency Q and energy (|12p are general and valid everywhere in the Galaxy. 
We have already numerically specified them for the region outside the bulge, where the velocity 
v o = V(r) = const. We will now discuss them for the central bulge. Note that only the frequency 
equation Q contains parameters r, JV* ~ M/ um /M , Mdm, vq that depend on the region in the 
Galaxy. The energy equation (fT2|) does not depend on the relative velocity of stars and PBHs vq 
because the actual velocity of collision is usually dominated by the escape velocity at the stellar 
surface. 

The frequency equation ^ scales with the relevant parameters as, 

ii ~ M lum M DM r~ 3 VQ 1 . (18) 

Let us now consider a model of the Galaxy in which the rotational velocity V(r) is proportional to 
the distance from the Galactic Center in the "bulge" and constant "outside" . This approximately 
follows what is observed. Distributions of the total (i.e. luminous plus dark matter) mass M to4 (r) 
and its density p(r) may be calculated from 

Yl * _ 4lGp = I^( r2 »), ^£ = 4 „v (19) 

r dr r z dr \ dr J dr 
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Results are given in Tabled One may assume that M; um (r) Mdm^) ~ M 2 ot (r). Then, in the 
bulge it is M 2 ot ~ r 6 , V(r) ~ r and therefore h ~ r 2 , which means that in the bulge the frequency 
of collisions is maximal for r = rr>] it does not increase towards the bulge center. Let us denote 
this maximal frequency by ho , and let us denote the frequency in the outside region by h out . From 
the above discussion and the last column of Table [2] it follows that, 

^ = i^4 = ^^ 20. (20) 

n ut M 2 out rj t r D 

In the view of our estimate at the end of Section 13.11 the collision frequency h ou t in the outside part 
of the Galaxy is extremely small. The factor of 20 does change much the situation in the bulge; 
the collision frequency hrj is also very small there. For this reason, in the rest of this paper we 
will not distinguish between the outside and the bulge, and take in our estimates of observational 
detectability n = n ou t as representative for the whole Galaxy. 

4. Signatures and detectability of a single encounter 

An encounter in which a PBH passes through a stellar object would be detectable if the following 
conditions are met: 

1. each single event is strong enough to be seen by the detector; 

2. the events are frequent enough so that they can be seen within the detector's life time. 

In subsection l4.lt we use the sensitivity of current detectors to calculate the maximum distance at 
which an event could be observed; in subsection 14.21 we estimate the required observation time for 
a detection. Both estimates are based on the energy loss and event rate calculated in Sections 13.21 
and 13.11 In Section [4.31 we determine the actual observation time required in order to see a single 
event in one of the current detectors. 



4.1. Detector sensitivity and photon fluence 

The aim of this section is to calculate the maximum distance at which an encounter involving a 
PBH passing through a Galactic object could be observed by a detector. This can be done by 
comparing the sensitivity of a detector, F sens , to the photon fluence from the event, F^. Detection 
is only possible for 

Fy ^ F sens . (21) 
The photon fluence and instrumental sensitivity can be determined as follows: 

• Photon fluence from a PBH passing through a Galactic object 

The photon fluence at Earth depends on the total energy deposited as given by Equ. (fl~3l) 



- 11 - 



and on the distance of the event from Earth, d, 

OSS /nn\ 

7 ~ n ■ (i + z) ■ <P ~ fi^ • 

The redshift factor (1+z) describes adiabatic energy losses, which can be neglected on Galactic 
scales, z ~ 0. The opening angle of the emitted signal Vl is determined by the shape of the 
shock. 

• Sensitivity 

For a positive detection, an instrument with an effective area A e ff needs to see at least one 
photon of energy EL, giving a fluence sensitivity of 

E 

Esens — ~~j • (23) 

This assumes that there is no additional background; as soon as there is some significant 
background, more than one photon would be needed for detection. We use this single photon 
argument to give an absolute lower limit on the fluence required in order for the signal to be 
detected. 



Combining Equations (|2Tj) . (f22l) and (|23|) . we have 



Eloss > E i (241 
fi • d? ~ A eff • 1 > 

Typically, the effective detection area A e ff depends on the energy of detection, and the detector 
is sensitive to the shape of the spectrum. Here, we optimistically assume that the effective area 
is concentrated at the lower energy threshold of the detector. For higher energies, the sensitivity 
would be less. For a given detector, the distance between the detector and the event to be observed 
must be less than 



Eloss 


A eff 


E'y 





is f^ » " ! *~ <25> 

Inserting into this the expression for the total energy released, as given by Equ. (|13j) . leads to 

d < d max = 4.7 pc • m pbhi20 ■ R 1 ^ ■ ^ — — j . p l Jl ■ A^ 250Q • n^ /2 ■ E~^ , (26) 

where values are given in terms of the X-ray detector XMM Newton: the effective area reference 
value is j4 e //,2500 = A e ff/ (2500 cm 2 ), at a reference energy of E 7j1 5 = E 7 /(1.5keV). The opening 
angle of the signal is expressed in terms of the maximum opening angle, corresponding to a signal 
emitted over half a hemisphere, 0,2tt = ri/(2 7r). 
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4.2. Observation time and interaction probability 

We can evaluate the observation time t a b s required for detecting the effect of a PBH passing through 
a Galactic object, by considering the event rate within the observable distance d max , n(d max ). The 
latter can be calculated from the total event rate, n, as calculated in Section I3.lt 

. / j x V (."max) . fnn\ 

n(d max ) = — n , (27) 

not 

where, V(d raa _ x ) is the volume containing the observable events 

V(d max ) =7rdl mx h (28) 

and 

V tot = 7T doalaxy h , (29) 

with (icaiaxy ~ 15 kpc being the radius of the Galaxy. We assume an isotropic distribution of sources 
within the Galactic disk with a height h = 300 pc. This means that we assume that the signal is 
always strong enough to reach Earth from 300 pc across the height of the disk. Realistically, some 
signals will be lost, since they are too weak to be detected at Earth, but as mentioned before, we 
work with the most optimistic assumptions in this paper. 

The number of events occurring within the observable distance d max from Earth is given as: 

N = n(d max ) • %^ • ^- ■ t obs . (30) 

4-7T 4-7T 

In this equation, the solid angle Q Q b s describes the Field of View (FoV) of the detector. The solid 
angle $7 accounts for a possibly beamed signal, which reduces the total number of events to be 
observed, since only a fraction Q/(4ir) is directed towards Earth. The total event rate n in the 
Galaxy was determined in Section 13. 1[ 

Using Equations ([27]) . ([28]) and (f29j) in Equ. ([30|) . the number of events within a radius (i max < 
^Galaxy is given as 

. /min(<i max ,<iG a iaxy)\ &>obs ^ . / olN 
N = n \ ) •l^I^-- (31) 

The radius of the Galaxy dcaiaxy gives an absolute limit for the number of events. For detection of 
a single event (N = 1), the required observation time is then 

4vr 4vr / n \ _1 /min (d max , Galaxy) \ ~ 2 , on , 
t obs = 1 yr • -— ■ — - - — — ■ 3 — . (32) 

il obs tt \lyr V V "Galaxy / 

Combining Equations ([9]), ([26]) and (|32|) gives 

U,K b ,<m^ k ) = 7.0-10 11 yr ■ m^ )20 ■ R~l ■ p~\ ■ ■ 

' ^ ' ^e//,2500 ' ^obs,0.2 ' E lA-5 (33) 
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for <i max < <i Galaxy and 

t ob s(m pbh > rr$*) = 1 • 10 3 yr • " ^ ' ^9 ■ T 1 

■ ^•^•^•( 12 -lOTcm/s ) (34) 
for <i max = (i Galaxy . Here, n o6S)0 . 2 = ^ o ;, s /(0.2 sr). 
This break in the function t b s (m p bh) occurs at a mass of 

break _o 9 , n 23 _ p-1/2 -1/2 (-,1/2 .-1/2 „l/2 /.m 

m pbh -3.2-10 g • • ^ — -1Q 7 cm/s J • P*,5 ' ^ • ^eZ/,2500 " #7,1-5 ' ( 35 ) 

The observation time reaches its minimum at this break mass. At lower masses up to the break mass, 
fnpbh < m pfeft, ak ' the observation time follows a rn~ bh — behavior. At higher masses, m p bh > m^l^t 
the observation time increases as fni bh . Observation conditions are therefore optimal at the break 
mass. 



4.3. Detectability 

As discussed in Section 13.21 the signals are expected to come at X-ray energies. We focus here 
on four detectors currently in operation which are well-suited for these investigations because of 
their relatively large FoVs and effective areas. First we consider XMM Newton wh i ch is s ensitive 
from extreme UV to X-rays. We use its effective area as quoted in lAschenbachl (|2002l ) at the 
reference energy 1.5 keV. The results which we obtain for this are equivalent to what is expected 
for observations with Chandra. For higher energies, from hard X-rays to soft gamma-rays, we use 
Swift-BAT and Fermi-GBM as reference detectors. Swift covers the range between 15 keV and 
150 keV; the high-energy detector of the GBM covers the energy range from 150 keV to 30 MeV. 
Both instruments provide quite large effective areas in combination with a large FoV, since they 
were designed to detect GRBs. For high-energy radiation, we use Fermi-LAT parameters, covering 
the range from 20 MeV up to 300 GeV. Table H] summarizes the properties of the four detectors. For 
Swift and Fermi, we use the lowe r energy thresholds as the reference e nergies for the effective areas 



quoted in (jBarthelmy et al.ll2005l : iMeegan et al 



2008 



Michelsonl 120031 ) . This is the most optimistic 



approach, since using higher energies would lead to reduced sensitivities. 

The detector's effective area, FoV and reference energy then determine the observation time and 
the break masses as discussed in the previous subsection. It turns out that the observation times 
for the four different detectors are related as 



t 



XMM 
obs 



Swift/BAT _ n Fcrmi/GBM _ „ „„„ Fcrmi/LAT 

1 ' Z obs — U - Z ' t obs — ■ t obs 



and the break masses m^ ak are related as 

break, XMM _ „ break, Swift/BAT _ „ . break, Fermi/GBM _ „ „„ break, Fermi/L AT 

m P bh — 1 ' m P bh — U.l • m pbh — U.U2 • m pbh 



(36) 



(37) 
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Swift/BAT gives the shortest observation time and is therefore best option concerning this. Fermi/LAT 
gives the largest break mass. Figure [2] shows observation time plotted against PBH mass for the 
case of XMM Newton in the relevant mass range of 10 15 g < m p bh < 10 26 g. Masses outside this 
range are already known not to contribute significantly (< 10%) to the dark matter. The solid line 
represents main sequence stars, the dashed line shows red giant cores, the dotted line is the result 
for white dwarfs and the dot-dashed line is for neutron stars. The shortest observation times are 
obtained for main sequence stars, and even there, the best case at m p bh = 2 • 10 g requires an 
observation time of more than 100 years. Similar results are obtained for Swift, Fermi-GBM and 
Fermi-LAT, as shown in Figures [21 [H and [SJ here as well, the observation times would need to be 
longer than 100 years in order to see one event, except in the case of Fermi/GBM at 3 • 10 25 g, 
where 60 years would be sufficient. The best observation times tH st are those for the break points 
in the plots of the observation time against PBH mass mj^. They are listed for the four detectors 
in Table H 

These observation times are far longer than the lifetimes of the present detectors and so, given 
that they are anyway optimistic lower limits on the real times, we conclude that we cannot obtain 
meaningful constraints in this way. The event rates are too low and the energy released is too 
small. 

Any constraint which we could give for the PBH contribution to dark matter in the Galaxy would 
be above 100%. We therefore do not convert the observation times into limits on the dark matter 
contribution. 



4.4. Alternative signatures 

There are a few other processes that could in principle lead to the production of a detectable signal 
for different detectors. We will briefly discuss here why these mechanisms also fail to produce a 
significant signal: 

• Shock acceleration of charged particles 

Particles can gain ener gy by scatte r ing o ff magnetic inhomogeneities via the mechanism of 



stochastic acceleration (|Fermilll949l . 1 19541 ). However, inside stars the densities are typically 
very high, and for white dwarfs and neutron stars the matter is degenerate. The mean 
free path of a proton in matter can be determined by considering the total cross section for 
proton-proton interactions (with a pp fa 50 mbarn): 

Am/p = — — — = 3 • 10" 4 cm • p* t5 . (38) 

P-k ' &pp 

This mean free path is typically smaller than a centimeter and does not allow for signifi- 
cant acceleration. Injection and acceleration of energetic charged particles in shock fronts is 
therefore highly disfavored, with the possible exception of the tenuous envelopes of red giant 
stars. 
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• Star quakes 

PBHs passing through an astrophysical object can induce star quakes. Although the emitted 
photons themselves cannot be distinguished from the general thermal spectrum of the star, the 
eve nt could be detected by Fourier-transforming the spectrum from time to frequency space, 
see (jDonea et al.lll999l ; iMartmez-Oliveros et al.ll2007l . e.g.) and references therein. However, 
these quakes can so far only be detected for the Sun and our rate for encounters between 



a PBH and the Sun is only 10 
excluded. 



yr . This method for getting detections can therefore be 



5. On the possibility of a PBH being captured by a star 



Roncadelli et al.l (120091 ) have considered the possibility of a PBH being captured by a star. In this 



section, we show that such a capture could happen only for PBHs with masses greater than about 
10 28 g, i.e. in the range already constrained by the micro-lensing results discussed earlier. The 
capture rates would therefore be reduced by a factor of 10 or more, considering that r/ < 0.1. 

The energy balance for a PBH being attracted by an astrophysical object in the formulation of a 
reduced two-body problem is given by 

E- m — E\ oss = E out . (39) 

Here, E\ oss is the energy loss undergone by the PBH as it passes through the star and E out is its 
final energy. Its initial energy is given by 

Ein = X • m pbh ■ Vq , (40) 

where we replace the reduced two-body mass by m^, because we are considering M* >> m p bh- 
The PBH can be captured by the star only if a sufficient part of the its initial kinetic energy can be 
dissipated ("lost") by the effective dynamical friction as it passes through the star. The condition 
for getting a bound system is E OVLt < 0. From Equations ([39]) and (|4"0]) we get 

Ei oss > E- m = - • m p bh ■ Vq . (41) 

This is the condition that needs to be fulfilled in order to have a capture. The energy loss from 
effective friction should be dominant, as determined from Equ. (|12p . It follows that 



2£i oss A _ 7 cm 2 ( v* V 2 „ (Equ.SU 



m pbh " 4 • 10 V ' m * 20 • U-lOWs ) " • P *' 5 V ° ■ (42) 

We can now derive a condition for the minimum PBH mass required for capture: 

mpbh > 20 > ^2.2.10°cm/s) ' (rW^p V ^ • ( 43 ) 



> ( r) ■ f 6.6 • 10 9 • • R~l + 1.2 • 10 7 ^ \ . r-i . p tU) 

V2.2-10 7 cm/sy V ' (2.2-10 7 cm/s) 2 y *' 9 ^ ' 
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The minimum mass for which capture can occur can now be computed for each of the astrophysical 
objects which we are considering (see table [3]): 

> 6 ' 6 - 1028 g (2.2. loL/^ +L2 - 1028g - (2.2 .toL/s)* ■»-» «!■—.(«) 
™* > ■ 1028 6 (2.2. ifaf - 1 - 2 - 10 " S (2.2. 10W.)* RG ° 8 ' <46> 

> 6 - 6 ' 10298 ' (2.2.1ot m /s^ +L2 - 102 ' g - ( 2 . 2 .lot m / S )' » hited ™ rfs - W 

2 4 

m„ fe/l > 6.6 • 10 27 g • ^ + 1.2 • 10 22 g ^ neutron stars . (48) 

P (2.2 • 10 7 cm/s) 2 (2.2 • 10 7 cm/s) 

On the view that the existence in the dark halo of PBHs with masses greater than about 10 26 g is 
already constrained by the micro-lensing observations, we conclude that the capture of a PBH by a 
star is difficult in any of the PBH mass ranges. Only relative velocities far below the average value 
could lead to a capture, but such events do not happen often enough to be significant. If any PBHs 
are captured by stars, this must be extremely rare, and certainly the fact that we do not observe 
them does not additionally constrain the PBH abundance in the dark halo of our Galaxj@. 

The above conclusion follows from the estimate of the energy lost when a PBH passes through a 
star, E\ OS8 , based on the Ruderman-Spiegel effective dynamical friction form ula. This may explain 



the difference between our conclusions and those of iRoncadelli et al.l (120091 ): they do not estimate 



£i oss , but instead assume that any PBH which collides with a star would be captured. 

We will examine the problem of estimating the probability of very low velocity PBH-star collisions 
at high PBH masses in a separate publication. 



6. Conclusions 

We have shown in the previous sections that the process of effective friction when a PBH passes 
through a Galactic object cannot lead to a significantly detectable signal with current instruments. 
We calculated lower limits for the observation time for the best-suited detectors: XMM Newton, 
Swift-BAT, and Fermi-GBM and Fermi-LAT. These lower limits were based on the assumptions 
that 

• The entire energy loss is radiated without significant delays and at the wavelength at which it 
is produced, 



2 A few years ago, Bohdan Paczynski suggested in a private conversation with M.A.A., that some unusual and 
very rare supernova-like events might be due to PBH collisions with white dwarfs. This suggestion was discussed 
later with Piero Madau, Rashid Sunyaev and others, but never followed up in detail. 
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• one photon is sufficient for having a significant detection of a signal, 

• all emitted energy is contained in the energy range of the detector, 

• the entire thickness of the Galactic disk can be observed, 

• the signal is emitted at the largest opening angle, i.e. 2ir. 

All of these assumptions were made in the sense of favoring detection; more realistic settings would 
reduce the detectability even more. We have investigated interactions with main sequence stars, 
red giants, white dwarfs, and neutron stars. In general, we can conclude that the best results are 
obtained for main sequence stars. However, independently of the class of objects, the difficulty 
is that either the total energy loss in the process is too small (for lower masses around 10 15 g < 
iTT-pbh < 10 20 g) or the event rate is too small (for larger masses around 10 20 g < m p f,h < 10 26 g). 
Even with the optimistic assumptions listed above, observation times of more than 100 years are 
required in order to see a single event for all of the objects and all of the detectors considered. 

Although we use average values in place of the mass and velocity distributions of PBHs, variations 
away from these central values will not influence the result of these calculat ions significantly. One 



could , for instance, argue that the clumping of dark matter as discussed by iBinnev and Tremaine 



(|2008l ) would lead to an enhanced signal. However, such an effect would need to increase the event 
rate by at least a factor of 10 2 in order to achieve reasonable observation times of t f, s < 1 yr. As 
shown in this paper, the central bulge of the Galaxy may improve the situation by a factor of 20, 
but this is still not enough to reduce observation times to realistic scales of less than a month. 

Hence, we conclude that constraining the abundance of PBHs in DM using interactions of PBHs 
with Galactic objects is not possible in the standard scenario described above. A significant rate 



of GRB-like events as suggested in (jZhilyaevi 120071 ) can therefore be excluded. There may be 



additional aspects, such as clumping of dark matter or observation of nearby dwarf galaxies, leading 
to an enhanced event rate and energy release. However, these will have to be able to improve the 
observation times in this standard scenario by at least a factor of 100 in order to give a significant 
signal. The main reason is that observation times of less than a month may be available with the 
different instruments, but not more. 

Concerning a possible collision with the Earth: the predicted rate is as low as 10~ 12 yr -1 • rn~ bh 2Q . 
This means that there would be less than one such event per 10 12— * 17 years for PBH masses above 
10 2 g. The interpret a tion o f the Tunguska event as a collision with a PBH of m v bh > 10 20 g, 
see I Jackson and Rvanl (|1973l ). can therefore be considered highly implausible on the basis of this 
straightforward statistical argument. The most common explanation discussed today is that a 
comet or asteroid disintegrated on its way to Earth. Since no macroscopic frag ments of the initia l 



object were found, its exact nature, comet or asteroid, is still unknown, see e.g. (jJopek et al.ll2008l ). 

On the other hand, PBHs with masses below m p f,h < 10 17 g may have struck or could still strike the 
Earth, as the encounter frequency may be greater than once during the Earth's life time. However, 
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the energy loss is too small to see such events by means of radiation. The possib ility of detecting 
them by means of an acoustic signal has been discussed bv lKhriplovich et al.l (|2008l ) but there it was 
pointed out that the signal is too weak to be seen by seismic detectors. It could be interesting to 
look for those signatures with the planned generation of high-energy neutrino telescopes, optimized 
to detect acoustic signals: large natural water, ice or salt reservoirs are to be i n strumented t o look 
for acoustic signals from Cherenkov radiation, see e.g. (jVandenbrouckd 120071 : iBeckerl l2008l ) for a 
review. The surface area of such a detector is planned to be of the order of 10 km 2 . If the PBH 
has a distinct signature, it might be observable with such an instrument. 
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Table 1. Basic properties of objects considered for collisions with primordial black holes. 
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Fig. 1. — Limits on the contribution to dark matter by PBHs. An absolute lower limit at 
around 10 15 g is given by the fact that PBHs with lower masses have evaporated by now. The 
limit at ~ 10 15 g is derived from the backgro und of photons at > 100 MeV in the Universe, 
see (jSreekumar et al.l Il998l ; iBarrau et al.l l2003l ) . Micro- lensing constraints in the mass range 
10 26 g< rupbh < 10 33 g, come from the EROS experiment (jAlcock et al.lll998l ). Limits in the mass 
range 10 33 g< rn p bh < 10 40 g, derived from measurements of the Cosmic Microwave Background by 



FIRAS and WMAP, are presented in lRicotti et al.l (|2008l ). In this figure, we only show the most 
stringent limits in the different mass regions and disregard those limits with weaker constraints. 



Table 2. Basic dependence of the velocity V(r), the mass density p(r) and the mass Mtot( r ) ° n 
the distance from the Galactic Center r. Here, Mtot{f) denotes the mass enclosed in the volume of 
the radius r. The parameters have different dependencies in the outer part of the Galaxy 

("outside") and the central part ("Bulge"). 
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XMM Newton (n„ ba = 0.2sr, A eff =2458cm 2 , E 7 =1.5keV) 




16 18 20 22 24 



l0 § [Mpbh/g] 

Fig. 2. — Observation time of PBH interactions with Galactic objects, observing with XMM New- 
ton. Solid line: main sequence stars; dashed line: red giant cores; dotted line: white dwarfs; 
dot-dashed line: neutron stars. The observation time decreases with m~^ h , until the point where 
the entire Galaxy can be observed, beyond which it starts increasing with m p bh- 




Fig. 3. — Observation time of PBH interactions with Galactic objects, observing with Swift-BAT. 
Same representation as in Fig. [2j 
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Fermi-GBM (n ob ,=27rsr, A.„= 1 2x 1 26cm 2 , E 7 =150keV) 
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Fig. 4. — Observation time of PBH interactions with Galactic objects, observing with Fermi-GBM. 
Same representation as in Fig. [2l 




Fig. 5. — Observation time of PBH interactions with Galactic objects, observing with Fermi-LAT. 
Same representation as in Fig. [2j 
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Table 3. Energy loss as calculated in Section [3T2l Equ. (fl"3]) for a PBH mass of m p bh = 10 20 jr. 
The velocity of the PBH at the surface of the star is calculated according to Equ. ([7]). The 
energy loss rate, dE\ oss /dt and the total energy loss £i oss , scale as (m p b/ l /10 20 g) 2 . The time it 
takes for the PBH to traverse the star is given by dt, see Section 1331 The peak energy E pea k is 
derived following Equ. (fT6|) in the case of regular stellar objects, and using Equ. (fT7|) for the 
relativistic objects (neutron stars and white dwarfs). The interaction probability n scales as 

("V>h/10 20 g) -1 according to Equ. ©. 



m pbh = 10 20 g 
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£ P eak [keV] 


0.2 


1 


11 
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h [yr -1 ] 


8 • 10 -13 
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11 


1.4 
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Table 4. Detector properties in the X-ray to soft gamma-ray range. In order to give an upper 
limit on the sensitivity, the lower energy threshold is used as E^ in the case of Swift and Fermi. 
For XMM Newton, the reference energy E^ = 1.5 keV was explicitly given for the listed effective 

area. 
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^min! -Etaax [keV] 

E 1 [keV] 
tt obs [sr] 
A eff [cm 2 ] 


(0.1; 15) 
1.5 
0.20 
2485 


(15; 150) 

15 

2 

5200 


(10; 3 • 10 4 ) 
150 
~ 2vr 
12 x 126 


(2 • 10 4 ; 3 • 10 8 ) 

2 • 10 4 

2 

8 • 10 3 


main seq. stars 
™" [g] 


2.1 ■ 10 24 
280 


0.7 • 10 24 
280 


2.1 • 10 25 
180 


1.4 • 10 26 
1500 


Red Giant Cores 
< [g] 

tot? [yr] 


2.1 • 10 24 
2.8 ■ 10 5 


0.7 • 10 24 
2.8 • 10 5 


2.1 • 10 25 

1.2 • 10 5 


1.4 • 10 26 

1.5 • 10 6 


White Dwarfs 

t b ot? [yr] 


4.2 • 10 24 
1.0 • 10 7 


2.1 • 10 24 
1.1 • 10 7 


4.2 • 10 25 
4.8 • 10 6 


2.8 • 10 26 
7.6 • 10 7 


Neutron Stars 

™S ak fe] 

t h ob f [yr] 


2.1 ■ 10 26 
2.8 ■ 10 4 


7.1 • 10 25 
2.8 ■ 10 4 


2.1 • 10 27 

1.2 • 10 4 


1.4 • 10 28 

1.5 • 10 7 



